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ABSTRACT 


The U.S. Marine Corps deploys Marine Air-Ground Task 
Forces (MAGTFs) by airlift or sealift to participactCemam 
numerous short-term exercises. These exercises are of such 
duration that resupply of the MAGTF by strategic airlift or 
sealift is not practical. Thus, only stocked spare parts 
are available for repairs during the exercise. 

A model is developed which provides the operational 
commander with a stockage policy for spare secondary 
reparables (e.g., tank engines, amtrack transmissions, etc. ) 
that optimizes the probability of successful mission 
completion subject to weight or volume constraints imposed 
by the MAGTF's mode of deployment. Optimization of ieee 
Stockage policy is stochastically modeled using data from 
the Marine Corps Integrated Maintenance Management System 
data base and then solved as an integer program. 

The integer program is coded using the Generalized 
Algebraic Modeling System language and solved using the 
zero/One Optimization Methods mixed integer program solver. 
Operational data for a Marine Amphibious Unit yields an 
integer program with 190 binary variables and 26 
constraints. A solution within 0.07% of optimality is 
obtained on an IBM 3033AP computer in 3.9 seconds and on a 


Zenith Z-248 personal computer in 176 seconds. 
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il NEROBUCTION 


The combat structure of the United States Marine Corps 
(USMC) centers around the concept of the Marine Air-Ground 
Task Force (MAGTF) which is a combined arms force consisting 
Seeoroungd, air, and combat service Support forces. The size 
of a MAGTF and its specific composition are dependent upon 
the particular mission assigned and the capability of the 
opposing forces. However, all MAGTFs are organized to take 
maximum advantage of the combat potential inherent in a 
rapidly deployable, closely integrated air-ground team under 
the control of a Single commander. A MAGTF consists of four 
basis elements: command element, ground combat element, 
aviation combat element and combat service support element. 
It is the mission of the combat service support element to 
provide all logistical support for the ground and aviation 
combat elements. [Ref. l:pp.1-4]. 

in this study, a model is developed which provides the 
Operational commander with a stockage policy for spare 
secondary reparables (rac, tank engines, amtrack 
Ghransmissions, etc.) that optimizes the probability of 
successful mission completion within the constraints imposed 
by the MAGTF's mode of deployment (i.e., airlift or 


sealift). The problem of determining this optimal stockage 


mp pies 
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policy is stochastically modeled using operational data 
resident in the Marine Corps Integrated Maintenance 
Management System (MIMMS) data base and then solved as a 
mathematical program. The remainder of this chapter will 
include a description of the scenarios which are of 
particular interest to the U.S. Marine Corps, an explanation 
of terms, and a definition of the problem. Additionally, 
justification for the approach taken here as compared to 


those in standard inventory models is provided. 


A. PROBLEM STATEMENT 

The U.S. Marine Corps currently deploys MAGTFs by 
airlift or sealift to participate in numerous short-term 
exercises such as BRIGHT STAR in Egypt, GALLANT EAGLE in the 
southwestern United States, and COLD WINTER in Norway. 
These MAGTFs range in size from 2,000 man Marine Amphibious 
Units (MAUS) to 16,000 man Marine Amphibious Brigades 
(MABS), depending upon the particular exercise. These 
exercises are of such duration that resupply of the MAGTF's 
ground combat element by strategic airlift or sealift from 
logistics bases in the continental United States is not 
preactical:. Furthermore, by doctrine [Ref. l:pp. 1-4], a 
MAGTF is task organized to maximize combat power; therefore, 
the number of mechanics, technicians, and other logistics 
personnel together with their test and repair equipment are 


kept to a minimum. Thus, the quantities of spare repair 


parts stocked by the combat service support elements are all 
that are available to support the requirements of the ground 
combat element for the duration of the exercise. 

In addition to these operational exercises, the Marine 
Corps deploys Marine Amphibious Units (MAUs) to the 
Mediterranean Sea as part of the U.S. Sixth Fleet ona 
continuous, six month rotating basis. These Mediterranean 
MAUS remain aboard ship except when conducting exercises. 
Therefore, the inventory of spares maintained by their 
combat service support element, called a MAU Service Support 
Group (MSSG), is limited by the cargo capacity of the 
assigned amphibious ships. Even though resupply from the 
continental United States by air is possible, it is desired 
to reduce dependency on this costly method free. csp. i]. 
Thus, Similar to the above exercise scenarios, the MSSG must 
stock sufficient spares, within the capacity pene tratnts of 
the assigned shipping, so as to satisfy the anticipated 
demands generated by the ground combat element's 
participation in exercises during the course of the six 
month deployment. 

Presently, a problem of even greater concern to Marine 
Corps combat service support organizations than either of 
those previously mentioned is the Marine Corps's role in the 
rapidly evolving maritime strategy of the U.S. Navy. As the 


amphibious power projection arm of the Navy, the Marine 
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Corps, in conjunction with the Navy, has developed and 
recently implemented the concept Known as the Maritime 
Prepositioning Force (MPF). According to doctrine [Ref. 3], 
a MPF consists of a Squadron of four to five ships loaded 
with combat equipment and supplies. One such MPF squadron 
is presently located in each of the following three ocean 
areas: Eastern Atlantic, Indian Ocean,and Western Pacific. 
When a MPF operation is ordered, a fly-in echelon consisting 
of 16,500 Marines and sailors, comprising a Marine 
Amphibious Brigade (MAB) and a Navy Support Element, is 
airlifted to a benign airfield in the vicinity of the 
objective area for link-up with equipment and supplies 
aboard the MPF squadron ships. The Brigade can be combat- 
capable and ready to move to an objective in five days or 
less and can operate for thirty days, independent of any 
Strategic airlift or sealift from the United States, uSing 
the supplies and Spare parts inventories on the ships [Ref. 
eo eee 6 Tete is this last requirement, to operate 
independently for thirty days uSing organic supplies and 
inventories, that poses the greatest challenge to Marine 
Corps combat Service support planners. However, as with 
Supporting short-term exercises and MAU deployments, the 
problem of supporting the MPF MAB simplifies to one of 
stocKing adequate spare parts inventories, subject to the 


capacity constraints of the MPF squadron ships, to meet the 


demand during the MAB's projected missions. tiie race chat 
such demand may vary unpredictably must be confronted. 

At the present time, to solve spare stockage problems 
Similar to those described above, Marine Corps combat 
service support planners determine the number of secondary 
items (referred to in Marine Corps' terminology as secondary 
Mepazables) to stock for each principal end item in the 
particular MAGTF's Table of Equipment (T/E) by uSing either 
last-period demand method or a simple moving average of past 
or historical demands generated during Similar exercises or 
Mediterranean deployments. In the situation where no 
historical demands are available, such as with the new MPF 
MABS, sStocKage levels are determined uSing estimated 
replacement rates in combination with the experience of 
maintenance and supply personnel. These techniques have 
proved inadequate [Ref. 2:p. 1] because they fail co 
consider such factors as changes to a MAGTF's Table of 
Equipment and variations in the operational schedule, both 
of which may occur from one exercise or deployment to the 
next. More importantly, they neglect to consider that the 
demands placed on the supply system are the direct result of 
equipment failures. 

To overcome these inadequacies, it is the objective of 
this study to develop a model which provides the operational 


commander at the MAU or MAB level with a stockage policy for 
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spare secondary reparables that optimizes the probability of 
successful mission completion within the constraints imposed 
by the mode of deployment. The proposed model is 
sufficiently general so as to be adaptable to any of the 
above scenarios and it can accommodate any additions or 
deletions to the MAGTF's Table of Equipment or variations 
in mission duration among any combination of principal end 
items. Additionally, the model incorporates a measure of 
effectiveness which provides the operational commander with 
an indicator of the overall impact of the stockage policy on 
Successful mission completion. Furthermore, the model 
utilizes the operational maintenance data resident in the 
Marine Corps Integrated Maintenance Management System 
(MIMMS) to determine operationally based failure rates. 
Thus, by use of this model a nearly optimal stockage policy 
for secondary reparables at the retail level can be 
determined for a wide variety of inventory problems 
currently facing Marine Corps combat service support 
planners. This will result in a more effective utilization 
of the limited space and weight resources available for 
spares’ inventories on shipping or aircraft assigned to 


deploying Marine Corps units. 


Be INVENTORY MODELS AND SOLUTION APPROACHES 

The stockage problems confronting Marine Corps planners 
can all be viewed as variations of the single period, 
constrained, multi-item inventory problem. The simplest 
such single-period model is the classic newsboy problem 
(Ref. 5], alternatively referred to as the Christmas tree 
problem. The characteristic distinguishing the models of 
this thesis from other standard inventory models, such as 
the classical economic-order-quantity or economic-order- 
interval models [Ref. 6:pp. 159-289] is that only a single 
time period, usually of finite length, is relevant; 
therefore, only a single initial procurement or order is 
made [Ref. 6:p. 297]. The previously described scenarios 
all fall within the context of sSingle-period models because 
the demand for secondary reparables is the direct result of 
equipment failures which occur at infrequent intervals 
during the course of the exercise. Resupply either does not 
occur or is kept to a minimum, and thus it can be ignored. 
Furthermore, it is Marine Corps policy to permit only 
functional unit replacement at the field or organizational 
levels of maintenance. Actual repair of secondary reparables 
at the intermediate maintenance level is very limited and 
for all practical purposes performed only at the depot 
level. Therefore, since the MAGTFs of concern in these 


scenarios are usually only authorized up to the intermediate 
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level of maintenance, the actual repair and return of 
secondary reparables to these retail levels of inventory is 
insignificant and can be ignored. This is particularly true 
when one considers that some of the exercise durations are 
of the same order of magnitude as the average repair time 
for some of the secondary reparables. Thus, each scenario 
consists of an exercise (or exercises) which is (are) of 
finite duration and during which repair and resupply are 
essentially nonexistent, so the combat service support 
planner must make a one-time decision of what to stock, not 
how much of each item to procure. He must do so within the 
constraints of the mode of deployment. This is analogous to 
the single period, constrained, multi-item inventory 
problem. 

There are several examples in the literature of single 
period inventory models applied to military logistics 
Situations. They include the spare parts kit problem [Ref. 
7:pp. 281-295], the fly-away kit problem [{Ref. 6:p. 324], 
and the submarine provisioning problem [Refs. 8:pp. 236-243, 
Oa tO lice The problem with each of these military 
applications is that they assume that there is some monetary 
stockout cost incurred when a demand exists and no spares 
are available to meet the demand. This is realistic in the 
Civilian sector where real monetary values can often be 


associated with each stockout. But, in most practaleaa 


feeeredry applications ft iS Gifficult, if not impossible, to 
determine a tangible, realistic, and objective stockout 
cost. This is particularly true for each of the previously 
described scenarios because there is no way to place a cost 
on the inability to complete a mission because of a 
deficiency of spares. ties, Lie rs GlttICult™ to apply” the 
models existing in the literature to these real-world Marine 
Corps inventory problems. 

This problem is not insurmountable. In the civilian 
sector when an organization does not know its stockout costs 
or is not confident in its estimates of them, management 
will set service levels which indicate the ability to meet 
customer demands from stock [Ref. 10:p. 149]. A service 
level is baSically the complement of the probability of 
stockout and really represents a subjective decision by 
management to accept some degree of stockout or customer 
disservice [Ref. 10:p. 150]. The U.S. Navy actually uses a 
Similar concept when establishing its Aviation Consolidated 
Aliowance Lists (AVCAL) and Coordinated Shipboard Allowance 
Lists (COSAL) [Ref. ll:pp. 4-30-4-50]. The service level 
measure is less than satisfactory for military applications 
Since it is subjectively established by combat service 
Support planners. service levels do not provide the 
Operational commander in the field with a meaningful 


indicator of the impact of the recommended stockage policy 


A xt 


REPRODUCED AT GOVERNNAIE il! F 


on successful mission completion. This is because there is 
no direct relationship between the probability of stockout 
and the probability of successfully completing a given 
mission. A high service level does not equate to an optimal 
probability of mission completion. This iS particu law 
true in the constrained, multi-item inventory situation 
where it might actually be better to accept a greater 
probability of stockout (lower service level) for a 
Secondary reparable with a low failure rate in order to be 
able to stock more of another secondary reparable with a 
high failure rate. 

To overcome the disadvantages associated with the 
minimization of stockout costs (actually variable costs) and 
service level approaches, the proposed model introduces a 


new measure of effectiveness, the probability of mission 


completion given the specified stockage policy, which eae 


be referred to as the probability of mission COmpiigiae 





throughout the remainder of thisS paper. This probability 
extends the familiar concept of reliability to include a 
system and its inva of spare parts [Ref. 8:p. 237]. A 
system is defined to be the quantity of a principal end item 
required to be operational throughout the duration of a 
given mission. This end item must be a “readiness 
reportable pacing item" as delineated in Marine Corps 


Bulletin 3000 [Ref 12: ppel=12 3a5e1 (Ze. A readiness 
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reportable pacing item is a principal end item whose 
operation is essential to the ability of a MAGTF to complete 
its assigned mission. Simply stated, the probability of 
mission completion is the probability that k-out-of-n (k<n) 
readiness reportable pacing items of the same type remain 
operational for a specified mission with a given stockage 
pPeLiLcy. Black and Proschan [Ref. 7:p. 283] proposed a 
Similar measure of effectiveness in their early work using 
marginal analysis to determine optimal spare parts kits at 
minimum cost. They called it probability of ecequaGu Or 
assurance of. adequacy, which is defined as 


The probability that either the equipment does not fail 


ewer toma. ™ -mabfunctiton or, rf it does fail for this 
reason, a replacement is available to provide continued 
successful equipment operation. [Ref. 8:p. 237] 


However, this definition applies only to series systems in 
which all equipment must remain operational for a successful 
mission. The proposed probability of mission wonptatron is 
more general since it can be extended to "k-out-of-n 
systems". This is required in the present study because 
the operational commander, not the military logistician, 
determines the quantities of a principal end item which he 
feels are necessary to complete the assigned mission. Ong 
example, a MAU commander may feel that he can accomplish a 
particular mission with just three out of five tanks which 
permits fewer tank secondary reparables to be stocked than 


if he required all five tanks to complete the same mission. 


-) 
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In this study the inventory problems in each of the 
previously described scenarios are modeled as single period, 
constrained, multi-item inventory problems in which the 
demands are the result of independent failures of secondary 
reparables. In the mathematical formulation of these 
models, the probability of mission completion is maximized, 
aS opposed to the standard approach of minimizing expected 
(variable) costs. This provides models which reflect the 
fact that in most military situations, and specifically in 
the three Marine Corps scenarios studied here, stockout or 
shortage costs cannot be determined. The inclusion and 
numbers of constraintS are dependent on the mode of 
deployment. However, in all cases studied, they are linear 
and will be referred to as capacity constraints. Thus, each 
model will be viewed as a mathematical program in which the 
objective function to be maximized is a nonlinear, 
nondecreasing function of the decision variable representing 
the quantity of secondary reparables to stock. Inclusion of 
the linear capacity constraint or constraints result in an 
overall integer nonlinear programming problem. Since the 
objective function is itself a complicated function of the 
decision variable, the solution approach is to linearize the 
objective function and then reformulate the entire problem 
using binary variables. This reformulation converts the 


original integer nonlinear programming problem into a binary 


IL 


linear programming problem which is then coded for computer 
processing using the Generalized Algebraic Modeling System 
(GAMS) [Refs. 13, 14, 15] and solved using the Zero/One 
Optimization Methods (ZOOM) [Ref. 16] mixed integer program 
solver. 

Since one of the primary purposes of this study is to 
provide the combat service support planners of the Second 
Force Service Support Group at Camp Lejeune, North Carolina 
with models solvable on a personal computer, it is desired 
to avoid dependency on commercial optimization packages not 
Mecally = Therefore, Ehe™~ initial solution 
methodology utilized Langrangian relaxation to establish an 
upper bound on the optimal solution of the integer linear 
problem. Then a heuristic based on marginal analysis is 
used to improve upon the best feasible solution obtained 
during the establishment of this upper bound, and this 
heuristically improved solution is accepted as the final 
solution provided that its associated objective value was 
reasonably close to the upper bound. This methodology is 
Similar to that implemented by DeWolfe [Ref. 17:pp. 28-32] 
to solve a Selective Reenlistment Bonus problem for the U.S. 
Marine Corps. The advantage of this approach is that it can 
be programmed in a language such as FORTRAN 77 thus avoiding 
dependency on a commercial optimizer. However, as the number 


of constraints increase the programming difficulty 


eS 


SPP UIGE 


REPRODUCED AT GOVERNMSIFII! F 


associated with the search procedure used in establishing 
the best upper bound increases, and for any more than two 
constraints actually becomes prohibitive. So, to keep the 
solution approach as flexible and general as possible, it 
was decided to solve all integer linear problems with the 
GAMS/ZOOM optimization package (commerically available in 
March 1988), thus allowing the solution of models with more 
than one constraint and avoiding restrictive programming 
considerations. 

The above solution approach overcomes the disadvantages 
associated with more standard solutions of the single period 
inventory model. For instance, it could be argued thac 
Since these inventory problems are all examples of decision 
making under risk, they could be solved using decision 
theory and payoff matrices to minimize expected costs. 
However, aS waS pointed out earlier, realistic monetary 
stockout costs are extremely difficult to determine in most 
military applications. Additionally, in the multiple item 
inventory problem, the payoff matrices become very large and 
the interrelationship among the payoffs associated with 
different combinations of items becomes difficult to 
evaluate and enumerate, so it is computationally more 
efficient to use a marginal analysis approach [Ref. 10:p. 
Zou Marginal analysis provides the basis for the 


techniques used in the spare parts kit problem [Ref. 
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ape zoOseZeemee end 1m the submarine provisioning problem 
(Ref. 8:p. 237], but it has practical limitations when more 
than one constraint must be considered. Hadley and Whitin 
[Ref. 6:pp 328-330] obtained exact integer solutions to 
their single constraint, fly-away kit problem using dynamic 
programming and Dreyfus and Law [Ref. 18:pp. 107-118] offer 
dynamic programming algorithms to solve single constraint, 
cargo-loading problems. However, they all point out that 
for problems involving more than two constraints even 
dynamic programming has computational limitations, 
Meticulariy when both the number of constraints and 
decision variables are large [Refs. 6:p. 331, 18:pp. 107- 
is |. The solution approach utilized in this study avoids 


these shortcomings. 


See lHESIS OUTLINE 
This thesis develops and presents a method for 
determining spare stockage levels of secondary reparables 


associated with the readiness reportable principal end items 


of a given Marine Air-Ground Task Force. In Chapter II, the 
maintenance data base, data analysis, and problems 
encountered during this analysis are discussed. The 


objective functions are defined, developed, and justified in 
Chapter III. Additionally, each scenario is modeled and 


formulated as an integer nonlinear programming problem for 
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which the objective function is linearized so that the 
problem can be reformulated as a integer linear problem to 
facilitate its solution. Specific solution methodology is 
presented in Chapter IV. The initial approach using 
Lagrangian relaxation with a heuristic is briefly described 
before the more general procedure using GAMS/ZOOM is 
developed and implemented. In Chapter V results uSing 
actual data for a typical Marine Amphibious Unit are 
presented and verification of the model is conducted using 
different sets of simulated failure data. The chapter 
concludes with a discussion of computational experience for 
both the Mainirame aneipersondal conpmeer. Chapter VI 
contains conclusions and recommendations. Finally, listings 
of both FORTRAN 77 and GAMS source codes are provided in the 


appendices. 


16 


Pree eo our AND ANALYSIS 


Peprimanry impetus for this study was the desire to 
develop a reliability-based, spares stockKage model for 
secondary reparables utilizing the maintenance data resident 
in the Marine Corps Integrated Maintenance Management System 
(MIMMS). Currently, spare stockage models within the Marine 
Corps are limited to those based on forecasts from past 
demands generated against the supply system. This 
forecasting approach neglects to consider that the source of 
the demands is the failure of equipment such as secondary 
-reparables, and it is precisely this type of failure 
information that can be extracted from relevant MIMMS sub- 
files which permits the development of a reliability-based 
model. 

AY DESCRIPTION OF DATA 

MIMMS iS an automated system which accumulates 
maintenance data on all serialized, principal end items in 
the Marine Corps inventory. It records and accumulates, in 
various sub-files, maintenance actions performed on a 
serialized end item during its service life at each of the 
three levels of maintenance (i.e., OEgGanrzatirtonal., 
intermediate, and depot) [Ref. 19]. In particular, the 


MIMMS sub-file of relevance to this study is the Equipment 
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Repair Order History File. This sub-file maintains all 
equipment repair orders (EROs) that were completed on every 
serialized end item during the past eighteen months, and it 
is updated on a quarterly basis. An ERO is opened each time 
a serialized end item is inducted into the maintenance cycle 
and contains information such as status of the end item 
(i.e., combat deadlined, safety deadlined, mission degraded, 
etc.), the meter reading at the time of the current 
maintenance action (expressed in either hours, miles, rounds 
or days depending Oe Ne end item), and the defect code 
which describes the nature of the failure. By looking at 
only those EROs in the MIMMS ERO History File which were 
opened for critical, combat deadlining repairs (1.e., 
category code M), a new, local user file can be created. 
This user file is organized by principal end item and 
includes only those end items which are considered 
“readiness reportable" [Ref. 12:pp. Encl (1) 1-23]. 
Associated with each end item are all serial numbers Og 
which a combat deadlining ERO exists; each serial number has 
a meter reading and defect code to reflect each ERO 
completed on that particular serial number during the 
eighteen month period covered by the ERO History File. 
Obviously, there may be multiple entries under the same 
serial number indicating that this end item suffered more 


than one combat deadlining failure during the reporting 
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period. There may be no entries if no failures were 
experienced by the particular serial number. Table 1 
provides a sample extract from this user-created data file 
and displays its contents for several tank serial numbers. 


TABLE 1 
EXTRACT FROM LOCAL USER DATA FILE 


METER DEFECT 

SERIAL # EOTC* READING CODE** 
502814 H 7a DOS (POWER TRAIN, CLUTCH) 
502828 H 51 A34 (ENGINE REPLACE) 

H 75 334 (COOLING SYS REPLACE) 

H 107 C34 (POWER RACK REPLACE) 

H 155 A34 (ENGINE REPLACE) 

H 163 B34 (TRANSMISSION REPLACE) 
502831 H 81 A34 (ENGINE REPLACE) 

H 149 B34 (TRANSMISSION REPLACE) 
502835 H 62 KO1l (ELEC SYS, GENERATOR) 


* EOTC= Equipment Operating Time Code 
(H=Hours, M=Miles, R=Rounds, D=Days) 


** Defect Codes are alphanumeric codes and have been 
annotated here for display purposes only. 
The data in this consolidated file can be used to estimate 
failure rates for each secondary reparable in either of two 
wayS and both methods will be described in the following 
section along with their advantages and disadvantages. 
It must be stated that both the format and method of update 
for many of the MIMMS files require improvement before a 
truly comprehensive analysis of the data is possible. A 
case in point is the ERO History File used in this study, 


and recommendations for its improvement are provided in the 
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following paragraph and in Chapter VI which minimize the 
impact on current MIMMS file management and the amount of 
effort required at the using unit input level. 

One of the basic problems with the ERO History File is 
that potentially useful data is being discarded during each 
Quarterly update in order to minimize file size. What is 
proposed here is to maintain a local (i.e., at the Force 
Service Support Group level) consolidated version of the ERO 
History File with a format similar to that described in the 
preceding paragraph. This file should be updated when its 
parent ERO History File is updated at the end of each 
quarter. Sinan, Chis consolidated file contains only 
information extracted from compa deadlining EROs, which 
comprise a small fraction of the total EROS in the parent 
file, there is a greater capability to accumulate relevant 
data without the need to automatically discard the oldest 
three months of data. The basic algorithm for each 
quarterly update is: 

Step 1. Select a combat deadlining ERO (Category Code M) 
from the ERO History File if one exists. If none 
exist, STOP. 

Step 2. Determine if the serial number associated with 
this ERO is already resident on the consolidated 
file. If it is, continue with Step 3. If ite 
not, add the serial number and its associated 
meter reading and defect code to the consolidated 


file under the appropriate principal end item. 
Return to Step l. 


20 


Step 3. Determine if the same defect code exists for this 
serial number. If it does, continue with Step 4. 
If it does not, add the meter reading and defect 
code to the list associated with this serial 
number. Return to Step 1. 

Step 4. Subtract the two meter readings to obtain a time 
between failure for this type of defect. Add 
this new data point to the list associated with 
this type of defect. Delete the record 
associated with the older meter reading (i.e., 
the smallest meter reading) and add the new meter 
reading and defect code to the end of the list 
associated with this serial number. Return to 
Step l. 

The result is an updated version of the consolidated file 
which is retained at the local level for update during the 
next quarterly cycle. Each readiness reportable principal 
end item in the file is broken down by serial number which 
in turn contain records of the most recent meter readings 
and defect codes. Each end item also consists of lists 
containing times between failures for each type of defect. 
Thus this local consolidated file does not alter any 
existing MIMMS file definitions, nor does it necessitate the 
creation of any new system file definitions. More 
importantly, it provides the combat service support planner 
With a consolidation of all the pertinent information 
required to compute failure rates and track the life history 
of any secondary reparable associated with a readiness 
reportable end item. This is a significant improvement over 


the present method which is restricted to using the most 


recent eighteen months of failure data resident in the ERO 
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History File, and it is this restriction tha” limics eo 
accumulation of failure times for certain types of secondary 


reparables (e.g., certain radio transmitters). 


By. DATA ANALYSIS 

Using the data resident in the consolidated file, the 
operating times between failures for each type of secondary 
reparable were computed. Failure rates were estimated using 
two different approaches. In the first approach, an 
exponential lifetime distribution was assumed for each 
secondary reparable and the maximum likelihood estimator 
(M.L.E.) for the parameter /; of an exponential distribution 
was used to provide an estimate of the failure rate. The 
second approach involves a BayeSian analysis of the times 
between failure to determine a posterior distribution for 
the failure rate of a secondary reparable. 


In deriving the M.L.E. for it is assumed that the 


Ae 
observed lifetimes of each secondary reparable i form a 
relevant random sample of size n; from an exponential 
distribution with unknown parameter ee The likelihood 
function must also account for the exposure times during 
which no failures of type i occurred. If these terms are 
not considered, the computed failure rate will tence. 


too large and excessive stockage levels would result. For 


the n; observed failure times and the observed exposure time 


a2 


feauring  =vneem no failures occurred, the likelihood 


mumctiOn 1S: 
ny 
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Maximizing the likelihood function with respect tov-, 


yields the following estimator of the failure rate: 


A MN; MN; 
A; — a —— . 
Ne ; 1h fee) 
ya i 
j=! 
where tT; represents the total exposure time for secondary 


reparable i. [Ref. 20:pp.282-296] 
The Bayesian approach also assumes that the distribution 
of the lifetimes for each secondary reparable i is 


exponential with parameter / however, it further assumes 


i’ 
that the exact value of this parameter is a realization of a 


random variable. Even though 4; is not precisely known, it 


has a prior distribution which is taken here to be a 
conjugate Gamma distribution with shape parameter 


Mere) and Scale parameter £.( / .>0): 


pri 
a =<) rama 
fe enon ie 0 


ee ep) = bla) (ae 
Q 0 . 





The prior distribution parameters are estimated using a 
parametric empirical Bayes (PEB) approach for failure rates 


as outlined by Gaver and Lehoczky [Ref.2l:pp. 220-224]. The 
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PEB approach uses the entire data set to compute estimates 
a ; and Ba which are substituted into the applicable 
formulas. To compute the prior parameters in PEB requires 


the maximization of the "marginal likelihood function": 


My 


La; B. Ne de Co) == a insti > 4imiim Bit 1a; =) —£ 3 * 
q 1) i a &) as | | ; Fel CXp T(a) Sim eXD é no - (233) 


i 


where, 


m individual source of data for component i (G@agme 
tank serial #502826) | 


\f; total number of data sources for component i 
failure rate for component 1 of data source m 


‘im observed exposure time of component i on data source 
Mm 


Nim Number of failures of component i on data source m 


oj shape parameter of Gamma prior distribution of 
failure rate for component 1 


b; scale parameter of Gamma prior distribution of 
failure rate for component i. 


Integration results in the following marginal likelihood 


function for secondary reparable 1: 


M, ‘ 
Lila;,B)3a;,6) = [ st ( ; “( = ™ (2.4) 
Ss ag LL tint Mad \ tim + By im +B; ) ‘ 


which is the product of individual negative binomial 


B 
distributions with parameters ee a and r=¢; : 
tia! ak 
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The numerical maximization of this function subject to Known 
upper and lower bounds on each «a; and f; is discussed in 
Chapter IV. The important result of this maximization is an 
estimate for the shape ( %;)and scale Glee parameters 
of the Gamma prior distribution that incorporates the actual 
observed operational data for each secondary reparable. 


Returning to the Bayesian analysis, the n; observed 


a 


, 


times between failure and the observed exposure times i 
Guring which no failures occurred are all a function of 


4; and give rise to the likelihood function, 


Ry 
. t+ ft’. : 
; Te _(%4 ‘ mae AT; eS 
ACs | Ai) = /; eXp j=l — 7 eX lee) 
where tT; once again represents the total exposure time for 
secondary reparable i. We are now interested in the 


updated distribution of 4, after the times between failure 


have been observed so as to provide a current estimate of 
future failure events. This is the posterior distribution 
and it is proportional to the product of the likelihood 


mumeerom (2.2) and the prior distribution (2.2): 


fA lt) c¢ fly lAadfp - 


Thus, 


f(A; 1 t) = cA exp 722 —! exp Bit 
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where c,; is a constant factor which depends on the observed 
times between failure, but does not depend on Aj The value 
of this constant can be determined uSing normalization. The 


resulting posterior distribution is: 


A A 
(c, 4+ py — _ en 
Aaj |) = fAj)data) = LE gad xp + By (2.6) 
[(n,; + x) 
a Gamma distribution with shape parameter n; +o; and scale 


parameter +B; . [Refv.s2Osap. 257=280) 

The posterior distribution could also have been derived 
directly using the fact that the family of Gamma 
distributions serves as a conjugate family of prior 
distributions for samples from an exponential distribution. 
Thus, using the theorem for Gamma conjugate priors [Ref. 
20:pp. 271-272], the posterior distribution for 4, Given 
the observed times between failure, is a Gamma distribution 
with shape parameter n; + «; and scale parameter 7+; . 
This agrees with the result derived in the previous 
paragraph and indicates that the posterior distribution can 
be defined knowing only the number of observations and the 
sum of all the exposure times. This iS an important result 
because it enables the posterior distribution to be updated 
very easily as new observations are obtained. 


Since the posterior distribution of A; has been shown 


to be a Gamma distribution, its expected value is: 
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This expected value also happens to be a "Bayes estimator" 


Em | data] —_ 


of «££: when a "Squared error loss function" is used [Ref. 


a 
BoD. 22/7). This estimator could also have been used to 
provide an estimate for each Age However, the primary 


purpose of the Bayesian analysis was not to provide a point 


estimate of + but rather to provide a distribution (i.e., 


eg 


posterior distribution) for each ee which could be used to 
make probabilistic predictions about demand functions 


mavolVvVing /.. This step is detailed in the next chapter. 
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III. MODEL FORMULATION ANDSEPES Che run 


In this chapter, the two different objective functions 
considered in this study are defined, derived, and 
justified. Additionally, each of the three scenarios 
described in the introduction is formulated as an integer 
nonlinear program whose objective function is subsequently 
linearized to create an integer linear problem. 

A. . MODEL DEVEEROCPRMENS 

The following formulation is developed to model the 
determination of secondary reparable inventory levels for 
each of the three scenarios described in the introduction. 
With the exception of the number of constraints, the 


formulation is the same for each scenario. 


INBECES 
i eer by secondary reparable or component 
DATA: 
Ww; unit shipping weight of secondary reparable i 
flbs] 
W weight capacity of aircraft or ship [lbs] 
Vy unit volume of secondary reparable [cuft] 
V volume capacity of aircraft or ship [cuft] 
uj maximum allowable number of spares for 


component i 
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FUNCTIONS : 


PA X;) probability that the number of readiness 
Bepomutable Oacing items (specified by the 
operational commander) containing secondary 
reparable i will survive a mission given x; 


aL 
Spares are stocked 
DECISION VARIABLE: 
2 Number of spares to stock for secondary 


reparable 1 


FORMULATION: 


x 
max PJs) Pr) 
i 


NV 
Sul, » wax, < IV 
a | 


i=!) 
\) 


Wye < 7 
a | 
i= 


OsnjSu i=1,..,N 


Seger 


The objective function is a mathematical expression of 
the probability that a quantity (specified by the 
operational commander) of each readiness reportable pacing 
item in a MAGTF's Table of Equipment will survive a mission 
of specified duration given a particular stockage policy for 
secondary reparables. This probability was referred to in 


the introduction as the probability of mission completion. 
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In terms of this definition, successful mission completion 
hinges on the survival of the MAGTF's pacing items. Each 
pacing item is composed of secondary reparables. The 
failure of any one secondary reparable will cause the pacing 
item to be combat deadlined and unable to complete its 
mission. Thus, each pacing item iS viewed as a system of 
secondary reparables in series. However, in a MAGTF's Table 
of Equipment, there is not just one of each type of pacing 
item, but rather a specific quantity of each type necessary 
to accomplish the assigned mission. This minimum quantity 
required will not only vary from one mission to another, but 
also among different MAGTF commanders when confronted with 
identical missions. eaaeee at the different pacing items 
taken as a whole do not constitute a simple series system in 
which all components must survive the given mission. 

In defining the objective function, two different 
functional forms were derived and studied. Both rely on the 
assumptions that secondary reparables fail at constant rates 
and that the failures of individual secondary reparables 
are independent of one another. The independence assumption 
is certainly true for the secondary reparables associated 
with different types of pacing items; however, it may not 
necessarily be true for the secondary reparables of a 
Specific pacing item. For example, the failure of M60Al1 


tank engines is independent of the failure of LVTP7AI1 
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amtrack engines, but the failure of amtrack final drives may 
not be independent of amtrack transmission failures. These 
assumptions together with the fact that secondary reparables 
from one type of pacing item are not interchangeable with 
those of another (e.g., a tank engine cannot be substituted 
for an eee engine) make it possible to consider only 
secondary reparables without having to specify the pacing 
items of which they are components. 

The functional form of the first objective function, in 
which the failure rate of each secondary reparable 1S 


assumed to be known, is as follows: 


qj Xj — Ki A 


, WD (ait _— 3.1) 
(x)= eee ex rod , 
PAX) ys an P 
a) 
where, 
a; number of pacing items containing secondary reparable 
a8 
k; minimum number of pacing items containing secondary 
reparable i required to be operational throughout 
given mission (k; < 4;) 
A 
AF M.L.E. of failure rate for secondary reparable i 
Fi mission duration during which secondary reparable i 


must operate. 
Mas function is a direct consequence of the constant 
failure rate assumption for each secondary reparable. 
Underlying this constant failure rate is an inherent 


exponential lifetime distribution which implies that the 
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number of failures defines a Poisson (arrival) process with 
exponential interarrival times. The number of failures of 
secondary reparable i to occur in any fixed interval of 
time, such as a mission of duration Be will have a Poisson 
e 


distribution with mean However, Since there are a 


Pas 
total of a; pacing items containing secondary reparable we 
and operating during this mission, the total number of 
failures will be the sum of a; independent, Poisson random 
variables (each with a mean aie) which is itself a Poisson 
random variable with a mean avt; . This is, really @ema 
an approximation since not all a; pacing items may be 
required during the Mission. In which case the overall 
failure rate would be smaller implying that we are probably 
overestimating the number of Spares required. On the other 
hand, it can be argued that with fewer than a, items the 
overall failure rate of the operational items may increase 
Since the same mission must be accomplished with fewer 
assets. In any case, this can be viewed aS a mild, worst 
case approximation. 

Therefore, equation (3.1) represents the probability 
that there are less than or equal to some specified number 
of failures of secondary reparable i during a mission of 
duration t;. This specified number is the sum of the number 


of spares, xX and the difference between the number of 


bd 


pacing items containing secondary reparable i (a;) and the 
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required minimum number of pacing items containing secondary 
reparable i (kK; ). As an example, consider the situation 
where the Table of Equipment of some MAGTF consists of five 
tanks each containing one engine (a;=5) and assume that the 
MAGTF commander decides that to accomplish the assigned 
mission he must have a minimum of four tanks operational at 
all times > Then, assuming one spare tank engine 
(X¥ =1) is stocked, there can be no more than two engine 
failures (a,+xX,-K,=5+1-4=2) to adhere to the commander's 
Guidance. It can be shown that this probability is equal to 
the probability that the quantity of each type of pacing 
item (or equivalently all its component secondary 
reparables) remains above a specified level (k,), because 
the following relation must always hold true: 
{number of failures of secondary reparable i} + 


{number of secondary reparables i operational and in 
inventory} = 


{total number of secondary reparables if, ee me 
Thus, 


P({number of secondary reparables i operational and 
in inventory at t}>k;) = 


P( {aj +X-} - {number of failures of secondary reparable 
a in dyox, i 


P({mumber of failures of secondary reparable 2 et } 
Saztx; IS hs 


Therefore, the probability of mission completion for the 


entire MAGTF is the product of the probabilities that the 
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individual secondary reparables exceed a level, k;, given x; 
Spares (or equivalently experience at most a pe 
failures). So, the maximum number of operational pacing 
items of a specific type at any given moment will be equal 
to the lowest level of its component secondary reparables. 
It is this minimum level that must exceed the specified 
level —- since if the minimum level exceeds kK;, all Temeie 
will exceed k;. Now, by the assumption that the secondary 
reparable failures are independent combined-.with the fact 
that if one type of pacing item cannot complete the mission 
the entire MAGTF will not be able to complete the mission, 
the probability that the minimum level exceeds the specified 


level (i.e., Dae) can be expressed as, 


a 
P) (x) } ¢ P>(x>) oe 8 Prlxy) = | [rte (om2) 


i=] 

which accounts for the functional form of the objeG@uawe 
functr ene 

Similar to the first objective function, the second 
objective function also takes advantage of the assumptions 
of constant failure rates for each secondary reparable and 
independence among the failures of different secondary 
reparables; however, it incorporates a more realistic 
approach as to the nature of the failure rates. In the 
previous objective function, the failure rate }),;, of each 


secondary reparable was assumed to be Known (actually an 
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estimate computed using the maximum likelihood estimator of 
the failure rate must realistically be used). For this 
second objective function, a Bayesian approach is taken in 


which the exact value of is explicitly recognized to be 


Ad 
unknown. The observed data (i.e., times to failure) are 
used to update the prior information concerning each failure 
rate to form a posterior distribution as discussed in the 
chapter on data analysis. This posterior distribution is 
necessary to derive the predictive distribution and the 
Pemetvonal tOrm of pl{x7), since we are interested in 
computing the probabilty distribution of the number of 
failures of secondary reparable 1 which occur during some 
mission of length t,;. By analogy to the previous derivation 
Sree (xs), 1t is desired to consider this number of 
failures because we want this quantity to be less than some 
eee teen fever (i-6.7 ay7*x--k,) sO aS to insure the entire 
MAGTF completes the assigned mission. Assuming that the 
Bibeare rate 47, is given, the total number of failures of 
Syeceenaary reparables of type i in a mission of length Cy 
will have a Poisson distribution with approximate mean 


re 


ay vq 1 


L 


; , F 2 (QjA jt; i =e “|e 
P(Exactly f; failures in (0,¢;]|/;) = exp _%*#4 


I} 


But 4; is not precisely known and has the following 


posterior density (as derived in the data analysis chapter) 


35 


Fale 4° 


FEE RP ROUUCED AT GOVERNKAIF ING F 2% 


given the observed times to failure and prior information: 


A A 
-+ Q: 
. (to) ty | 
f(i;| data) = ——- 2 FT! exp 4ilti t BD 
A 
ae) 
where, 


ty total exposure time for secondary reparable i 


n; number of observations for secondary reparable 


x 

oF empirical Bayes estimate of shape parameter for Gamma 
(odes Cr Dee eDELOG OLS tai eon 

A 

6, empirical Bayes estimate of scale parameter for Gamma 


( J,.0;,8; ) prior distribution. 


This can be rewritten as, 


mit Biles By agp Akt + Bi 
IV | data) = (A(t + B))* " \(t,+ B)) exp A(t; + Bj) 
i ae eee 
C(n; + «;) 


l 
Now, to calculate the predictive distribution, 


P(Exactly f; failures in (0,1;]) =| P(Exactly ff failures in (0,1) 14) flz;\datalds; 
0 


which is upon substitution, 


: : - esl ae A 
co, (aA ,t,V! } ime) ear: 1 eres | 
| (SEE exprayty( SRP pig l exp AT PND, . 
0 Ji: D(2; + «;) 


Evaluation of the integral yields: 
A A A 

P(Exacty:f failures in (0,0])=——L eee | rrr )" oS ( ___ a } _ 

fl U(n; + &,) at +t, + Bi ate 7, ae 

This is a negative binomial distribution with parameters 


A 
tT; + B; A 
2. a) ae and =n +; . 
ait; Seis or B; 
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Even though p,;(x;) is defined as the probability that 
the stockage of x; spares for secondary reparable i allows 
the pacing item of which it is a component to exceed some 
minimum level ( en); it has been shown that this 1s 
equivalent to the probability that the number of failures of 
secondary reparable i during mission t,; is less than or 
Pa@iial tO .a,+x,;-k; ~ (i.e, P(number of failures of secondary 


reparable i in (O,G5/< de) i Therefore, the second 


BumetLTOnal form of pets) is: 


a; + x; — k; A a an oe 
D(n; + a; + m) Tie pe Dep Ce Qt; (3,3) 
pls) a - A ; Zn ) 
ay et ee) ajt;+ 1, + B; at; +t, + B; 
where, 
a; number of pacing items containing secondary reparable 
‘i 


k; minimum number of pacing items containing secondary 
reparable i required to be operational throughout 
given mission (k; <a; ) 


Nn,  mnumber of observations (data points) of times between 
failure for secondary reparable i 

mi total observed exposure time for secondary reparable i 

; mission duration during which secondary reparable i 


must operate 


a. empirical Bayes estimate of shape parameter for Gamma 
(oencmSp de PotOnedistmi bution 


empirical Bayes estimate of scale parameter for Gamma 
Oe) ) Prior distribution . 


The reasoning behind the functional form of equation (3.2) 


BevVomewee 101 the first pr(x,) Stzil holds true for the 
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newly derived p;(x,;); therefore, the functional form of this 
objective function is still the product of the individual 
pi@t,) s- 

The linear weight capacity restricts the total shipping 
weight (i.e., secondary reparable plus container) for all 
secondary reparables stocked to be less than or equal to the 
weight capacity of the ship or airplane. Similarly, the 
linear constraint on volume or space capacity requires that 
all secondary reparables stocked not exceed the volume 
available. The final constraint in Pl indicates that there 
is an upper bound on the number of spares for each secondary 
reparable that can be stocked. This iS a very real 
constraint because there are only a finite number of each 
type of secondary reparable at the wholesale inventory level 
from which the retail inventories used to support all 
exercises, deployments, and contingencies can be 
established. Thus, a combat service support planner must 
consider the requirements of the force as a whole so as not 
to commit too many secondary reparables to any one 
exercise. 

The number of constraints is dependent on the scenario. 
Any scenario that involves the airlift of the ground combat 
element and its supporting forces into the exercise or 
Operation area, such as the short-term exercises described 


in the introduction, will have only one constraint. This is 
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the constraint on the allowable cargo load (ACL). Allowable 
cargo load is: 


Pre ins cral linewting taccor in determining airlift 


requirements. .. . This is a weight limiting factor 
expressed in terms of STON [short tons] which cannot be 
exceeded in planning aircraft loads. The ACL is 


established by the Air Force for each type aircraft for a 
specific mission. [Ref. 22:pp. 4-9] 


Since the formulation of problem Pl now reduces to just one 
constraint, it is recognizable as a nonlinear knapsack 
problem due to the nonlinear objective function. The 
individual secondary reparables represent the items or 
Pemmodities that must fit into a knapsack while the ACL 
constraint corresponds to the kKnapsack's weight or volume 
Capacity. 

The other two scenarios described in the introduction 
both involve two constraints. These two constraints 
represent the weight and volume capacities of the shipping 
assigned to perform the sealift. In the case of the MAUs 
deploying to the Mediterranean Sea, the secondary reparables 
are normally stored aboard one of the ships of the 
amphibious task force such as a LHA, LPH, LKA, LSD, or LST. 
Thus, the constraints reflect the weight and volume capacity 
of the ship's cargo hold or compartment which has been 
designated for secondary reparable storage. On the other 
hand, in considering secondary reparable support of the MPF 
MAB, the secondary reparables are spread-loaded among 


several of the MPF squadron ships; therefore, the weight and 
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volume constraints are actually the combined weight and 
volume capacities of the cargo holds of the specified MPF 
squadron ships. Since both of these sealift scenarios 
involve two contraints, problem Pl can no longer be regarded 
aS a nonlinear knapsack problem, but it is still an integer 
nonlinear programming problem. 

Regardless of the number of constraints or whether 
equation (3.1) or (3.3) is used as the functional form for 
Da (ee the objective function is a complex, nonlinear 
function of the decision variable. For this reason, problem 
Pl will be transformed into a binary linear programming 
problem. The following section explains this transformation 


and the new formulation. 


B. TRANSFORMATION TO A BINARY LINEAR PROGRAM 

It can be seen from the final constraint of problem Pl 
that the decision variable, 2 eee es constrained to a set of 
discrete values. Therefore the objective function will 
first be separated using the properties of logarithms and 
then the transformation technique for converting discrete 
variables to binary variables will be used to reformulate 


the problem *[ Ref 2 Sele | = Taking the natural logarithm 


of the objective function yields: 


N NV 
In (Tren) 7 >In Pp Ax 4) ° 
i=] i=l 
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Now, the final constraint in Pl can be expressed as, 


x; = >; = (Si, Sits ele. Sigs} = One oe u;} 


and this is equivalent to the constraint set [Ref. 23:p. 


ez): 
Us 
a 

Xp PSP 

i=0 


i; 


é 


oe el 


J=U 
Ny oa my — Oe uy 2 


Applying the separation technique and transformation to 


problem Pl results in the new formulation: 


PNDICES ; 
=|, ee secondary reparable or component 
| es number of secondary reparables 1 
DATA: 
i in(p;(j)) 
ay total shipping weight for j} spares of type i 
[lbs] 
W weight capacity of airplane or ship [lbs] 
“i total volume of j spare of typei[cuft] 
V volume capacity of airplane or ship [cuft] 
ay maximum allowable number of spares for 


component i 
DECISION VARIABLE: 
ae { if j spares of type i are selected 
| 0 


otherwise 
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FORMULATION: 
N UW; 


max » » SipXij (ee) 
t=1 ;=0 


AN i 
sy: ; 
st.) ) wy SH 


aii] S {Q, 1} 


Assuming Known failure rates, the coefficients, Laz: in the 


objective function are: 


wei Kk; A ad) A 
v (a;/.t;) pe 
fe=inf > exp (3.4) 
‘ ass AL. 
m=0 
where, 
a; number of pacing items containing secondary reparable 
iL 
kK; minimum number of pacing items containing secondary 


reparable i required to be operational throughout 
given mission (kK, <a, ) 


A 
4; M.L.E. of failure rate for secondary reparable i 
fi; mission duration during which secondary reparable i 
must operate. 
If perfect knowledge about the failure rates is not assumed 


and the Bayesian approach is adopted the coefficients are: 
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tas _ A 
“Sn + 9; + Tn; + oj + mt) Gree one ( ail; (3 
faoed ( - ) : ) 


Aa ‘ 
m=0 m: ae si a) Qt; ni T; cu B; Ql; ote i a B; 


where 

Q; number of pacing items containing secondary reparable 
ai: 

ki; minimum number of pacing items containing secondary 
reparable i required to be operational throughout 
given mission (kK; <a; ) 

Nn; number of observations (data points) of times between 
failure for secondary reparable 1 

t; total observed exposure time for secondary reparable 1 

iF mission duration during which secondary reparable i 
must operate 

2; empirical Bayes estimate of shape parameter for Gamma 


eee b) Pre rmemotrd bution 
f; empirical Bayes estimate of scale parameter for Gamma 
ate ae Creel St Er bu tT On. 
It is of interest to note that for situations in which 
the airlift model is applicable, problem P2 is just a 
special case of the generalized assignment problem. In the 
formulation of the generalized assignment problem 
[Ref. 24:pp.345-346], each agent (or man) can perform or be 
assigned more than one task (or job) provided that the 
resource available to the agent is not exceeded; whereas, in 
the classical assignment problem [Ref. 23:pp. 61-62], each 
agent can be assigned to only one task. IGE |sienele 


formulations, each task must be assigned to exactly one of 
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the agents. The airlift model can be interpreted as a 
special case of this generalized assignment problem in which 
one agent must perform all N tasks, however, he has (uy; + 1) 
alternative ways to accomplish each task. 

Problem P2 is equivalent to Pl and even though the 
number of variables and constraints have increased, it is 
the preferred formulation because it eliminates the complex 
nonlinear objective function of Pl. The coefficients ij 
can be directly computed using recursive formulas for the 
Poisson or negative binomial cumulative distribution 
functions; whereas, p,(x,;) could not since it was itself a 
function of the decision variable. Additionally, because of 
the assumption of the independence of failures, other 
failure distributions can be readily introduced into this 
linear formulation. Furthermore, the formulation is general 
enough so that additional constraints can be incorporated as 
required and as will be seen in the following chapter, it 
lends itself easily to programming in an algebraic modeling 


language. 
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IV. SOLUTION METHODOLOGY AND IMPLEMENTATION 


The purpose of this study is to determine an optimal 
stockage level of spare secondary reparables for each of the 
three different scenarios described in the introduction. 
Of course "optimality" is with reference to the failure 
model adopted and the relevance of the data used to estimate 
its parameters. The determination of the stockage levels 
is the subject of this chapter. The binary linear nature of 
the formulation presented in the preceding chapter suggests 
that several different methods of solution are available to 
determine the optimal stockKage policy for a given scenario. 
Each alternative method will be discussed and justification 
for the method chosen iS provided in addition to its 


implementation. 


A. SOLUTION APPROACHES AND METHODOLOGY 

The mode of deployment in a particular scenario 
determines the number of constraints. The initial solution 
approach to the scenario in which the airlift model was 
applicable utilized the technique of Lagrangian relaxation. 
Although it is standard procedure to establish bounds for 
integer programs (such as problem P2) uSing a linear 
programming (LP) relaxation, it was initially desired to 


avoid dependency on commercial LP solvers not locally 
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available to the combat service support planners of the 
Second Force Service Support Group at Camp Lejeune, North 
Carolina. Therefore, Lagrangian relaxation was adopted 
because it lends itself nicely to programming in readily 
available FORTRAN 77. It must be emphasized that the 
structure of problem P2 did not demand that this specific 
solution technique be applied. In fact, the LP relaxation 
of problem P2 could have been solved using an LP solver. 
Thus, Lagrangian relaxation was utilized not out of 
necessity, but primarily because it offered a realistic 
alternative to the use of a commercial LP optimizer. 

Briefly stated, the technique of Lagrangian relaxation, 
as Outlined by Fisher [Ref. 25:pp. 10-21], involves moving 
complicating constraints into the objective function using 
the product of the Lagrangian multiplier and the constraint 
violation as a penalty term. In particular, examination of 
problem P2 in the context of the airlift scenario results in 
the following binary linear earl ee (actually a 
generalized assignment formulation): 


INDIGES: 
i=1,..,N secondary reparables or components 


j=9,.-,4; number of secondary reparables i 
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DATA: 


Sij in(p;(j)) 

a total shipping weight for j spares of type 1 
[1bs} 

MG allowable cargo load of airplane [1bs] 

u; maximum allowable number of spares for 


component 1 


PeclLsSLlON VARIABLE: 


room fe ostscomen=cype 1 are selected 
= U0 otherwise 
FORMULATION: 
NY 
ae NON + (P3) 


2,2, SRG 
i=| s=0 


N UU; 
Vay 
Au Hg VA aah < af Ge 
f=] 7=0 
u; 
yxy a Bel 
| 
j=9 
xij = {0,1} 


iaesenly complicating constraint is the allowable cargo load 
waxes) constraint. The relaxed formulation with scalar 
Lagrangian multiplier, 4 , consists of N separable, multiple 
choice problems which are easy to solve. On the other hand, 
for those situations in which both the weight and volume 
constraints of problem P2 must be considered (i.e., 
Mediterranean MAU or MPF MAB scenarios), the Lagrangian 
multiplier, 4, is now vector-valued (i.e., 4=( 41, 45) since 


two constraints must be relaxed, and this is no longer a 
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trivial problem to solve. The primary reason for this 
complication is that the solution approach now involves a 
two-dimensional search which is more difficult to program 
than the simple bisection search associated with the single, 
scalar Lagrangian multiplier. The problem is not impossible 
to solve and methods exist, such as a coordinate search 
method suggested by DeWolfe, Stevens, and Wood [{Ref.26: 
pp.8-10], which appears to offer good results for this type 
of problem. 

As additional constraints (such as budget, holding 
(storage) cost, or shortage (stockout) cost constraints) are 
added, the programming difficulty associated with the search 
procedure of the Lagrangian function becomes prohibitive. 
Therefore, to accommodate the possible inclusion of 
additional constraints in the formulation of problem P2, it 
was decided to utilize the Generalized Algebraic Modeling 
System (GAMS) in conjunction with the Zero/One Optimization 
Methods (ZOOM) mixed integer program solver. The use of 
this optimization package not only avoids the restrictive 
programming considerations discussed above, but it can also 
incorporate the PEB as an initial step, prior to solution of 
PEP Oormre . 

The approach to the PEB is to maximize the marginal 
likelihood function associated with each secondary reparable 


1 with respect to a; and §$,;. This marginal likelihood 
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variables. The —LoOmrmulatton tor ChiS MaximiZation problem 
is: 
rie LCES: 
i=1,..,N secondary reparable or component 
Pil.) sourcesot Gata (e.g., tank serial number 508282) 
DATA: 
imp number of observations of component i failures 


for data source m 


im total observed exposure time for component i of 
data source m. 


Oj upper bound on the shape parameter of the Gamma 
Paton GdastribuivaoneLor COMNpOmenia 1. 


B; upper bound on the scale parameter of the Gamma 
Pater ALStrebution for component i. 
DECISION VARIABLES: 


Oj shape parameter for Gamma prior distribution 
associated with compenent i. 


B; scale parameter for Gamma prior distribution 
associated with component i. 


FORMULATION: 
MM; A A A 
aie at) B; =i lim "im 
hee) | (Pa) 
Nim" I (a.;) eats Bb; ee p; 
m=! 


Wee 0<a;< 4G; 
A — 
0< Bp. <B; 
Standard maximization by taking partial derivatives with 


respect to xk; and f; does not yield a satisfactory closed 
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form solution which can be practically implemented. Thus, 
this nonlinear maximization problem is also solved using 
GAMS in conjunction with the nonlinear program Solvers 
MINOS [Ref. 27]. These values are then incorporated into 
the optimization of problem P2. In summary, under the 
current solution approach, the entire model can be coden for 
computer processing on either a mainframe system or personal 
computer and solved using a Single modeling system (i.e., 
GAMS) and its compatible solvers. This approach is 
Sufficiently flexible to incorporate additional constraints 


Or decision variables. 


Be IMPLEMENTATION 

Since the decision was made to determine the optimal 
stockage policy for all scenarios uSing GAMS, no specifics 
of the original solution approach involving Lagrangian 
relaxation will be mentioned other than the fact that it is 
implemented by means of a FORTRAN 77 computer program. The 
FORTRAN program is based upon a Similar one developed by 
DeWolfe [Ref. 17:pp. 23-24]. A copy of the source code 
developed for this study is provided in Appendix A and 
examples of its two required input files and of the single 
Output file are given in Appendices B and C, respectively. 

The solution approach adopted centers around the coding 
of problems P2 and P4 for computer processing using the GAMS 


language. The reader is referred to the excellent tutorial 
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by Rosenthal [Ref. 28] for the specifics of this versatile 
mathematical modeling language. Upon coding, the problems 
are solved using the GAMS compatible solvers of MINOS and 
ZOOM for nonlinear programs and mixed integer programs 
respectively. Specifically, the source code provided in 
Appendix D incorporates the successive solve feature of GAMS 
to first solve the nonlinear maximization problem 
represented by problem P4. The results of this maximization 
provide the estimates of the shape (a. ) and scale (Bs) 
parameters of the prior distributions which are required 
before the second binary linear maximization (P2) can be 
performed. The final result is the optimal stockage level 
of spare secondary reparables for the given scenario. Even 
though the source code in Appendix D is specific to the 
airlift model (P3), GAMS makes it very easy to incorporate 
additional constraints. For example, if a volume constraint 
is also applicable, a single GAMS PARAMETER statement to 
accomodate the volume data and EQUATION statement to reflect 
the volume constraint are the only additions necessary. 

The source code provided is compatible with either 
mainframe or personal computer usage. When run on the IBM 
3033AP mainframe computer, GAMS version 2.05 was utilized in 
conjunction with MINOS Version 5.1 and ZOOM version 2.1 and 
on the Zenith Z-248 personal computer (equipped with 80287-8 


math coprocessor and 640 Kbyte memory expansion board), PC 
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GAMS Version 2.04 was utilized in conjunction wath” tie. 


Version 5.0 and ZOOM Version 2.1. 


a 


Ve. WRESULTS 


This chapter summarizes results obtained using actual 
Operational maintenance data to determine an optimal level 
of spare secondary reparables for a Marine Amphibious Unit 
(MAU) in a realistic exercise scenario. Verification of the 
model is accomplished uSing Simulated failure data. 
Estimation of the parameters of the prior distribution by 
an analysis of the data is discussed and demonstrated; this 
realizes computationally the PEB procedure. Emad 7, 
computational experience for both the mainframe and Zenith- 
248 personal computer is briefly discussed. 
ae MODEL TESTING 

The data used in the development and testing of all 
models was provided by the Operations Section of the 
Supported Activities Supply System (SASSY) Management Unit 
in coordination with the Maintenance Information Systems 
Wemagetmeme Jrtice Of the Second Force Service Support Group 
at Camp LeJeune, North Carolina. This data was extracted 
from the ERO History File of the Marine Corps Integrated 
Maintenance Management System and analyzed according to the 
procedure presented in Chapter II. This analysis provided 


the basic parameters for equation (3.5) (or (3.4)) from 
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which the objective function coefficients were generated 
using GAMS. 

The airlift model (P3) was used to determine the optimal 
stockage policy for twenty-five different secondary 
reparables associated with the pacing items of a Marine 
Amphibious Unit participating in a realistic short-term 
exercise. This problem generates 190 binary variables and 
26° CONStLalnes. The recommended stockage level for the 
secondary reparables associated with the three pacing items 
(i.e., M60A1l Tank, LTP7A1l Amtrack, and M1O1A1 105mm 
Howitzer) in this exercise is presented in Table 2. Table 2 
displays the levels computed using both the GAMS/ MINOS/ZOOM 
package and the Lagrangian relaxation technique. 

For this exercise scenario 99.89% of the allowable cargo 
load was utilized and the objective function value was 
within 0.07% of the upper bound established during the 
branch and bound procedure used by the ZOOM solver. The 
probability of mission completion given the stockage level 
was 0.4149 and this highlights one of the advantages 
associated with using this model. This low probability 
implies that the allowable cargo load must be increased. The 
Only way to accomplish this ites to increase the number of 
cargo planes assigned to this exercise. This provides 


Marine Corps planners with a tangible basis on which to 
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TABLE 2 
MAU AIRLIFT SCENARIO SPARE STOCKAGE LEVELS 


SPARE STOCKAGE LEVEL 


GAMS / LAGRANGIAN 
NSN END ITEM ZOOM RELAXATION 
2815001245387 M60A1 Ze Z 
2520002241867 M60A1 S S 
mel5002395819 M60A1 Z Z 
Zo20001549632 M60A1 2 Zz 
Zoo0000886657 M60A1 2 S 
2520001184942 M60A1 2 Z 
Zao lO13580 2 M60A1 il iD 
merOO1TO0703803 M60A1 il i 
mO25019388105 M60A1 4 5 
ieebOO10720397 M60A1 il uF 
meZ2010708993 M60A1 nH ih 
2815004303480 CVE AE S 4 
mm29003973384 LVTP/7A1 5 > 
2520001443385 LVTP7A1 4 4 
5805014593410 LVTP7A1 4 2) 
2520008949532 LVTP7A1 Z Dp 
Zoe0O011714636 LVTP7A1 2 S: 
mao0011509757 Evrrvae Z eZ 
5865012207848 LVTP7A1 S S 
2320002035660 LVTP7A1 Z Z 
220002317276 iy LP 7 Ae He i 
2530000886650 LVTP7A1 Zi 7b 
mega LL855059 M101A1 5 4 
1005011086434 M1IO1A1 = 5 
1005011457709 M1O1A1 5 S 
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justify their requests for additional airlift assets fren 
the Military Airlift Command of the U.S. Air Force. Table 3 


Summarizes the results for both solution approaches. 


TABLE 3 


SUMMARY OF RESULTS FOR MAU AIRLIFT SCENARIO 


Objective Probability Deviation 
SOLuEron Funct2en of Mission Percentage from 
Technique Value Completion of ACL Used Optimaliigy 
GAMS /MINOS / 
ZOOM =O eo, 0.4149 99.89% 0.07% 
Lagrangian 
Relaxation -Q0.8820 Ov4 iss 100% 0.47% 


(w/Heuristic ) 


B. “~MODEB  VEREPECAT ION 

Verification of the modeling approach was conducted by 
using the model to verify the fact that the negative 
binomial predictive function reduces to the Poisson-Gamma 
function as the number of failures observed (and accordingly 
the total exposure time) approaches infinity. Thus, the 
Poisson function provides the limiting case or best possible 
stockage level attainable assuming that the failure rates 
for each secondary reparable are precisely Known. 

To accomplish this verification, the MAU airlift model 
waS run uSing the Poisson objective function (3.4) with 
fixed failure rates. The resulting stockage levels 


represent the best stockage policy that can be achieved. 
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Then, using those failure rates, twenty-five different sets 
of simulated failure data were generated for each of 10, 25, 
50, 100, 200 and 500 failures. These simulated failure 
times were randomly generated using the probability integral 
transform to obtain exponential failure times from uniform 
(O,1) deviates provided by the LLRANDOMII Random Number 
Generation Package [Ref. 29]. The Bayesian analysis of this 
data was performed resulting in the negative binomial 
predictive distribution and the objective function, equation 
ere) ) . A noninformative Gamma prior distribution hice. | gaa’ 
and [ =1 was employed in this process. StocKkage levels for 
this method are displayed in Table 4. Note that about 100 
failure observations are needed before the stockKage levels 
approach those computed using the actual failure rates 
(i.e., the optimal stockage policy given perfect 
information). Figure 1 is a plot of the interquartile ange 
(IOR) of the probabiltiy of mission completion associated 
with each failure sample size. This graph demonstrates that 
as the number of observed failures increases, the 
variability (represented by the IQR) in the probability of 
mission completion decreases. 

In addition to the above verification, the fact that the 
model and solution approach responded predictably to a wide 
range of different failure observations and exposure times 


indicates that the model and approach tend to be robust. 
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MAU AIRLIFT SCENARIO SPARE STOCKAGE LEVELS 


NSN 


2815001245387 
Zac0002241367 
ZelsOUZS7oG 13 
2520001549632 
2530000886657 
2520001184942 
ZIZOOTOES S302 
1010010703803 
102501 2333205 
LOTOCRS7 Zs 97 
10250107089¢3 
2815004303480 
2520003973384 
2520001443385 
9805014593410 
2520008949532 
2910011714636 
Z53 0011505727 
9865012207848 
4320012035660 
Z9Z0900Z31 7275 
2530000886650 
1005011855059 
1005011086434 
1005011457709 
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TABLE 4 


SPARS STOCKAGE Gly oe 


NO 
) 


OOO PRrRPRNNNNNOhRAPOrFRORFPRPERFNN A W 


NEGATIVE BINOMIAL 


50 


OUWRRNNNNNOPROARORFP RRP RFPRPRNN bh W 


98 


LOS 


OU PRrERNWNNNOLAARRPRF RPRPRNNWNWN 


200 


OQUOPRFRNWNNNOhVUIPRPRP Rr RRPRNN WNW NH 


900 


OU PNRrFNWNNNOROARRPRPRPRPRPRNWN WN 


POISS@GRH 


OO PNRFNWANNNORORRPRrPPRPRPRPNWNWNDN 


O.60 


COS 


0.45 


PROBABILITY OF MISSION COMPLETION 
0.40 


0.go 


! ] ! | _ 


100 200 300 400 Seo 
NUMBER OF FAIRURES 


Figure 1. Variability in P(Mission Completion) with 
Number of Observed Failures 
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GC. BAYESIAN AND EMPIRICAL BAYES ANALYSIS 

Bayesian analysis and particularly PEB suggest 
themselves when the failure rates of different units (e.g., 
tank engines on different tanks or amtrack transmissions on 
different amtracks) are sufficiently alike, but individual 
experience is small, so that "borrowing strength" by a 
Suitable pooling process can improve individual estimates. 
In this study, the failure rate associated with each 


secondary reparable ( ) was assumed to follow a Gamma 


Ad 
prior distribution. This prior distribution represents the 
probability that the unknown value of the failure rate lies 
in various regions of the parameter space [Ref. 20:p. 260]. 
In other words, incorporating the prior distribution enables 
a tighter bound on the unknown value of the failure rate to 
be obtained than would be possible by using the data 
pertaining to that rate alone. The primary disadvantage of 
the ordinary or classical Bayesian approach is that the 


prior distribution's parameters (i.e., a. and 6 are 


ae i? 
assumed to be known. This deficiency is overcome in this 
study by employing a PEB approach, in which the entire data 
set ( (fj. lin) M=l,..,3M; ) associated with each secondary 
reparable i is used to estimate the parameters of the prior 
distribution. These estimates are then used in the standard 
Bayesian analysis to compute the posterior and predictive 


distribucrens: The results are unavoidably influenced by 
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the particular choice of prior used (here a Gamma). 
Sentitivity tests have not been conducted, but indication is 
that the procedure is rather insensitive. 

A specific indication of this insensitivity is seen by 
comparing the airlift model stockage levels obtained when 
minimal prior information is included to the levels obtained 
when the empirical Bayes estimates of the prior distribution 
are included. The results of this comparison are summarized 
in Table 5. 

The levels in the column labeled "NONINFORMATIVE PRIOR" 
were obtained by running the airlift model with the randomly 
generated failure times (described in the previous section) 
associated with the number of failures sample size of ten. 
A "diffuse" Or noninformative joe al oe G@iastribution 
(i.e., %4;= [;=1) was used. A diffuse prior distribution is 
"informationless in a relative sense and means only that it 
is diffuse relative to the sample information". [Ref. 30:p. 
198] The same airlift model was then run using empirical 
Bayes estimates for the prior parameters and these levels 
appear in the column labeled "INFORMATIVE PRIOR". These 


levels are only slightly closer to the Poisson levels’ than 
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NSN 


Z2E8US O01 245387, 
ZIZOUOZZALSS/ 
2815002395819 
2320001549532 
2530000865657 
2920001184942 
ZIZOCRO RS SZ 
1010010703803 
LOZS029S sees 
POROORO 2 OS97 
10Z50T07039S2 
2815004303480 
2920003973384 
2520001443385 
9805014593410 
2920008949532 
2910011714636 
ZIS00 TS O97 a7 
5865012207848 
4320012035660 
ZIZQOOZS1 7276 
2930000886650 
MGO5011855059 
1005011086434 
1005011457709 


TABLES 


EMPIRICAL BAYES ESTIMATION 


STOCKAGE LEVELS 


NEGATIVE BINOMIAL 


NONINFORMATIVE 


PIenok 


> G) 
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POISSON 
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are those computed when minimal prior information is 
considered. 

Even though these stockage levels are only slightly 
closer to the Poisson levels for this specific instance, by 
simply using the available data for 10 failures to compute 
empirical Bayes estimates of the prior distribution 
parameters, spare stockage levels comparable to those 
obtained in the previous section using failure sample sizes 
of between 50 to 100 can be achieved. This is particularly 
Significant since we have no control over the number of 
failures observed for each secondary reparable or their 
exposure times and would not uSually observe this number of 
failures (i.e., 50-100) during the reporting period covered 
byes the ERO |Histony File. Also, since operational 
reliability data may not be available for the exact location 
where the MAGTF is deploying, pooling the data from sources 
which have operated in similar environments has the 
advantage of providing more realistic failures rates than if 


this prior information were not considered. 


D. COMPUTATIONAL EXPERIENCE 

So as to Keep the implementation and solution of the 
models developed in this study as general and flexible as 
possible existing software packages (i.e., STATGRAPHICS, 
GAMS, etc.) were utilized. The goal was to obtain solutions 


in a relatively user-friendly personal computer environment 
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requiring a minimum of operator interaction. heretic. 
computational efficiency in terms of computer run time (or 
algorithmic complexity) was not of primary concern and is 
included here solely for completeness. No attempt was made 
to employ special algorithms that exploit problem structure. 
This is not to downplay the importance of this factor, but 
rather to emphasize that during this study the primary 
concern was to obtain solutions to real-world problems as 
easily and efficiently as possible with a minimum of 
programming effort. 

Running time for the airlift model on the IBM 3033AP 
using the GAMS/MINOS/ZOOM package was 3.9 seconds in 
contrast to the 176 seconds obtained when run on the Zenith 
Z-248 PC. It is possible on the IBM 3033AP to combine the 
PEB calculations with the maximization of problem P3 by 
uSing the successive solve feature of GAMS; however, this 
cannot be done on the Zenith Z-248 PC since internal 
parameters in the PC version of GAMS are exceeded. 
Therefore, each of the empirical Bayes maximizations must be 
performed separately when using this PC. The resulting 
estimates have to be manually input into the appropriate 
GAMS PARAMETER statements before maximization of problem P3 
can be accomplished. 

Even though the Lagrangian relaxation procedure with 


heuristic was not adopted as the solution technique for the 
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reasons discussed in Chapter IV, it provides much faster 
Selution times to the airlift model (P3) on both the IBM 
3033AP and Zenith Z-248 PC. The Lagrangian procedure, coded 
in FORTRAN 77, provides a running time of 0.32 seconds on 
the IBM 3033AP and about 30 seconds on the Zenith Z-248 PC 
where it was coded using Microsoft FORTRAN (Version 3.2). 
Finally, comparison of the stockage levels obtained 
using GAMS/MINOS/ZOOM with those obtained using LP 
relaxation for each of the 150 Simulated failure time cases 
revealed that in every case the LP relaxation provided 
identical stockage levels with the exception that it 
fractionated for exactly one of the twenty-five secondary 
reparables. So, if a general policy was adopted to round 
down to the lower of the fractionated stockage levels (So as 
not to violate the ACL contrast), the LP relaxation would 
provide comparable results to those obtained using the mixed 
integer program solver, ZOOM. This 1S an important result 
because the run times associated with the LP relaxation were 
Significantly less than those using ZOOM. ineparcticular, 
if the rounding down policy is applied to the LP relaxation 
of the airlift model in this study, the resulting integer 
stockage level provides an objective function value which is 
within 1.87% of the LP upper bound or 1.73% of the upper 
bound established by the branch and bound procedure utilized 


by ZOOM. This is not quite as good as the result obtained 
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Girectly from ZOOM (within 0.07% of optimality), however 
for those situations in which computer run time may be a 
consideration (e.g., as problem size increases) or in which 
very tight bounds may not be required, a very good solution 
can be obtained from the much simpler and less costly LP 


relaxation. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


In this thesis, the problem of determining spare 


stockage levels iO Ie deployable MAGTFs was 
solved by 
(a) stochastically modeling demand utilizing operational 


data from the Marine Corps Integrated Maintenance 
Management System, and then 


(op) using mathematical programming to optimize a 
Suitable, operationally relevant objective function. 


Memaccomplisch this, a new measure Ofmeitfectiveness was 
defined which provides the operational commander on the 
ground with an objective indicator of the overall impact of 
the stockage policy on the probability of mission 
completion. This measure of effectiveness was evaluated as 
the end result of a probabilistic modeling step, a 
statistical data analysis, and a subsequent optimization. 
The generality of the model formulation makes it 
adaptable to any of the three scenarios of present interest 
to the Marine Corps. The use of a PC-compatible algebraic 
modeling language (GAMS) makes it particularly easy to add 
Or delete constraintS as required. Specifically, the model 
and solution approach were applied to a realistic exercise 
scenario and an essentially optimal stockKage policy was 


determined which maximized the probability of mission 


Oye 
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completion within the constraints imposed by the deployment 
mode. 

During the analysis of the MIMMS data base (1.e.,ERO 
History File) numerous problems and inconsistencies were 
encountered. As a result, the estimates of the failure 
rates for some of the secondary reparables are higher than 
would be the case if all the potential times to failure 
could be extracted from the data base. To overcome this 
discrepancy, it is highly recommended that the quarterly 
update procedure described in Chapter II be conducted at the 
Force Service Support Group level. If this cannot be done, 
then at the very least the meter readings associated with 
all pacing items should be recorded and input into MIMMS 
prior to each quarterly update. This would insure that 
beginning and ending meter readings are available for each 
serialized pacing item during the reporting period covered 
by the ERO History File and this would prevent the loss of 
valuable exposure time data. 

The limited strategic lift assets available to deploying 
Marine Corps units demands that effective utilization be 
made of these resources. Since only a small percentage of 
these assets can be allocated for spare parts inventories, 
emphasis must be placed on stocking those parts which would 
have the greatest impact on the completion of the assigned 


mission. The reliability-based model developed in this 
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study accomplishes this objective by maximizing the 
probability of mission completion within the constraints of 


the mode of deployment. 
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APPENDIX A 


LISTING OF SOURCE CODE FOR LANGRANGIAN RELAXATION 


PROGRAM RELAX2 
CHARACTER NSNX¥13,NOMENX7 


INTEGER NCOMP,NEQUIP,NSPARE,MAXSPR, IDNUM, BSPARE, SPBND,NDENS, 
XNOPER,NXPOSE 


REAL X8 P,W,MAXWT,DINF,OPTDEV, ZLOW, ZUP,PSUC,HILAMD,WI,MTBF, 
XDURMIS, ZHLOW, XPOSE,ALPHA, BETA 


PARAMETER CNCOMP=25,NEQUIP=3,MAXSPR=10,MAXWT=44400.0D0, 
*DINF=1.0D20) 


DIMENSION NSNCONCOMP),NOMENCNCOMP) ,NSPARECNCOMP),IDNUMCNCOMP), 
X¥PCHCOMP,0:MAXSPR),WCNCOMP,0:MAXSPR), BSPARECHICOMP),SPBNDCNCOMP) , 
XHDENSCNCOMP) ,NOPERCHCOMP) ,XPOSECNCOMP) ,NXPOSECNCOMP), 
XDURMISCNCOMP),ALPHACNCOMP), BETACNCOMP ) 


CALL INIT2C(DINF,MAXSPR, NEQUIP,MAXWT,P,W,NSN, NOMEN, RHILAMD, NSPARE, 
XSPBND,NCOMP, IDNUM,NDENS,NOPER, XPOSE,NXPOSE, DURMIS,ALPHA, BETA) 


CALL BOUNDCDINF,P,W,NCOMP,MAXSPR, SPBND,MAXWT,HILAMD,NSPARE, ZLONW, 
XZUP,VWT, BSPARE)D 


CALL HERISTCDINF,P,VI,NCOMP,MAXSPR, ZHLOW,WT,MAXWT, SPBND, BSPARE, 
XLARRAY) 


WRITE CLOGSaD 
50 FORMAT ('l]'/'0!,!---------------------------------------------- 7 
XP ee ee ee t) 


OPTDEV = -1.0D0 X* 100.0D0 * (1.0D0 —- (Zneonaz)) 


WRITE €10,100) ZHLODI 
100 FORMAT (1X, "OBJECTIVE FUNCTIONS VALUE Gi Se Deemer 
WRiE Cae0, 1950) OFMDEW ; 
150 FORMAT (€'0',."'THIS SOLUTION IS WITHIN *,D17.10,' % OF THE OP isi 
PRINT *, 'ZUP=',ZUP, 'ZLOW=',ZLON, 'ZHLOW=', ZHLOW 
PSUCT="EArCChoo 
WiLL C10 27 5) see 


175 FORMAT ('O', "PCMESSLON SUCCESS) = "aye 
WRITE €10,180) 

180 FORMAT ('0!, !---------- 33-2 nn on 1 
0 ae ee y 


WRITE €10,200) 
200 FORMAT €'0',5X,'NSN*',10X,'END ITEM',6X, "SPARE LEVEL") 


DO 250 I=1,NCOMP 
WRITE €10,300) NSNCI)D,NOMENCIDNUMCI)),BSPARECI) 
250. CONTLEUE 
300 FORMAT €'0',A13,7X,A7,7X,19) 
PRINT %, "WEIGHT USED=',WT, 'MAXWT=',MAXWT 


Sic 
END 


70 


SUBROUTINE INIT2CDINF,MAXSPR,NEQUIP,MAXWT,P.W,NSN,NOMEN,HILAMD, 
*¥NSPARE,SPBND,NCOMP, IDNUM,NDENS,NOPER,XPOSE,NXPCSE, DURMIS,ALPHA, 
XBETA)D 


INTEGER NCOMP 
REAL 8 MAXHWT 


INTEGER IDNUMCNCOMP),NDENSCNCOMP),NOPERCNCOMP),LMT,MAXSPR, 
XNSPARECNCOMP),SPBNDCHCOMP),NEQUIP,NXPOSECNCOMP ) 


REAL XS WONCOMP,0:MAXSPR),DURMISCNCOMP), DINF,P1,Q,BETACNCOMP), 
eben CHeoMr, O:MAXSrR),ATLAMD,HLTEMP,CUMSUM, XPOSECNCOMP ),PP, 
XALPHACNCOMP ) 


CHARACTER NSNX13CNCOMP) ,NOMENX7 CNCOMP ) 


Rene Gi P00 SCHRSNCI I DNUMCT),XPOSECID,NXPOSECI),WC1,1),SPBNDCI), 
XALPHACI)D,BETACI), I=1,NCOMP) 
19 FORMAT (A13,2X,12,2X,F9.1,2X,19,2X,F6.1,2X,13,2X,F/.4,24%,FE.2) 


READ (2,20) CNOMENCI),NDENSCI),NOPERCI)D,DURMISCI)D, I=l1,NEQUIP) 
20 FORMAT (CA7,49X,13,4X,13,9%,F6.1) 


DG cl vl =13ascOMP 
PRINT X,NSNCI),IDNUMCI),XPOSECID,NXPOSECT)I,NOI,1),SPBHNDCI), 
XALPHACI)D, BETACTI) 
at CONT INUVE 


DO 22 T=1,NEQUIP 
fo ons nomen CLI,NDENSC I; NOPER CI, DURMISCI) 
ec CONTINUE 


Doe 1 =i, NCOMP 
NSPARECI) = 0 
DO 50 J=0,MAASPR 
PCIl,J) = DINF 
eS GNSte J eNiGl, J) = Pili 
50 CONTINUE 


25 COM THUE 


DS 75 T=1,NCOMP 
Pl = (XPOSECTI+BETACIJI/SCNDENSCIDNUMC I) DX 
AUUMGmoOCLONUNC lL) © APOSEC IT) + BETACI)) 
Cy 000 = Pl 
pipes ee CHXPOSECI) + ALPHACTI)) 
Pigeons CDDNUMCII) = NOPERCIDNUM( 13) 
DO 100 J=0,SPBNDiI) 
IF (J .NE. 1) WCI,J) = JRWCTI,1) 
CUMSUM = 9.0D0 
BOP e> K=-O,CEMT + J) 
EP Cpe o: 0) THEN 
Pai. = 1.000 


ZbS (2 
PieMe | eR am ct 2 TKO TIAL. OD0I7RIZG 
END I 
Cwiisvre—- CUMSUM + PTEMPI 
125 CONTINUE 
Cursum = CUMSUM * PP 
mes? = EOGCECUMSUM) 


rn ee, b= py 073, EN=*.PC1,3), 'CUMSUM=',CUMSUM 
100 COMPIHUE 
75 CONTINUE 


Gill 
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ea> 
200 


50 


HILAMD = -DINE 
DOe 200, T=1 NCOP 
DO 225 J=); sre 


HLTEMP = ( PCIS2) - Pll, 0) Jo” sie 
IF CHLTEMP .GT. HIGAMDO MEAD = ee 
CONTINUE 
CONTINUE 
PRIN? *, "HILAMBDA=!,HILAMD 
RETURN 
END 


SUBROUTINE BOUNDCDINF,P,W,NCOMP,MAXSPR, SPBND,MAXWT,HILAMD, NSPARE, 


XZLOWN, ZUP,WT, BSPARE)D 


INTEGER NCOMP,BSPARECNCOMP),MAXSPR,SPBNDOCNCOMP) ,NSPARECNCOMP ) 
REAL 8 BSTLAM, ZLOW,WT,POCNCOMP,0:MAXSPR) .WONCOMP,0:MAXSPR),DINF, 


XMAXWT,HILAMD, EPS,LFTEND,RTEND,ZLBEST,LAMBDA, ZUP 


Reo = HEEAMDZT 00046 


_LFTEND = 0.90D0 
RTEND = 1.1DOXHILAMD 
ZLBEST = -DINE 
LAMBDA = CLFETEND + RTEND)/2.0D0 


CALL MAXFCTCDINF,P,W,NCOMP,MAXSPR, SPBND,MAXWT, LAMBDA,NSPARE, ZLOW, 


*¥ZUP,WT) 


IF (WT .LE. MAXWT) THEN 
RTEND = LAMBDA 
IF C4e0hee GE. ZUBESl) THEN 


ZLBEST = ZLOW 
BSTLAM = LAM3DA 
pie ir 
jae sal 
LFTEND = LAMBDA 
Eel F 


IF C CRTEND = LRTEND) (261 Er Sacer ta 
CALL MAXFEACDINF,P,W,NCOMP,MAXSPR,SPBND, BSTLAM, BSPARE,WT, ZLOW) 


RETURN 
END 
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SUBROUTINE MAXFCTCDINF,P,W,NCOMP,MAXSPR,SPBND,MAXWT,LAMBDA,NSPARE, 
XZLOW, ZUP,WT) 


INTEGER NCOMP,NSPARECNCOMP),MAXSPR, SPBNDCNCOMP), INDEX 


REALX8 LAMBDA,ZLOW,WT,PCNCOMP,0:MAXSPR),WCNCOMP,0:MAXSPR),DINF, 
XCTOTAL,CMPMAX, ZVALJ,OBJFCT,OBJVAL,ZUP,MAXWT 


WT = 0.0D0 
CTOTAL = Q0.0D0 
ZLOW = C.QODO 


DO 25 I=1,NCOMP 
CMPMAX = -DINF 
INDEX = 0 


DO 50 J=0,SPBND(1) 
ZVALJ = PCI,J) - LAMBDAXWCI, J) 
OBJVAL = PCI,J) 
IF (ZVALJ .GE. CMPMAX) THEN 
Nemec )) = J 


CMPMAX = ZVALJ 
OBJFCT = OBJVAL 
INDEX = J 
END .iis 
50 CONTINUE 


CTOTAL = CTOTAL + CMPMAX 
PU pemeeOnet OBJECT. 
NT = WT + WCI,INDEX) 

a, CONTINUE 


ZUP = CTOTAL + LAMBDAXMAXWT 
RETURN 


END 
Sa MAK CROCUS pp NeCUME Aono Sho, Sol LAM, BSPARE,WT, 
¥ZLOV) 


INTEGER NCOMP, BSPARECHCOMP),MAXSPR, SPBNDOCNCOMP), INDEX 


REAL 8 BSTLAM, ZLOW,WT,PCNCOMP, 0:MAXSPR),WONCOMP,0:MAXSPR),DINF, 
XCMPMAX, ZVALJ,OBJFCT, OBJVAL 


Wr = 0.0D0 
ZLOW = 0.0D0 


DO 25 IT=1,NCOMP 
CMPMAX = -DINF 
INDEX = 0 


DO 50 J=0,S5PBNDCI) 
ZVALJ = PCI,J) - BSTLAMXWCI, J) 
OBJVAL = PCI,J) 
IF (ZVALJ .GE. CMPMAX) THEN 
BSPARECI) = J 


CMPMAX = ZVALJ 
OBJFCT = OBJVAL 
INDEX = J 
Eid) 1 F. 
50 CONTINUE 


ZLOW = ZLOW + OBJFCT 
Wl = WT + WCI,INDEX) 
25 CONTINUE 


RETURN 
END 


fee 
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£2 


50 


fie, 


100 


125 


SUBROUTINE HERISTCDINF,P,W,NCOMP,MAXSPR, ZHLOW,WI,MAXWT,SPBND, 
XBSPARE, LARRAY ) 


INTEGER NCOMP, BSPARECNCOMP ),MAXSPR,SPBNDCNCOMP), INDEX 


REAL X8 DINF,PCNCOMP,0:MAXSPR),WONCOMP,0:MAXSPR),WT1L,WT,MAXWT, 
XXSSWT, CHKVAL,LARRAYCNCOMP), NUMER, DENOM, LBEST, ZHLOW 


CHKVAL = DINF/1.1D0 
MTL = WT 
XSSWT = MAXWNT - WT] 


DO 50° f=1,NCOMnP 
IF CCBSPARECI) .GE. SPBNDCI)) .OR. CPCl1, GBSPARECT) +1) ieee 
X¥CHKVAL)) THEN 


LARRAYCI) = -DINF 
GOTO 50 
ay Bee a 
NUMER = PCI, ¢CBSPARECI)+1)) = RCI, BSPeke 
DENOM = WCI, CBSPARECI)D +1009 =" 2 Soe ae Oe 
IF (€CDENOM .GT. XSSWT) (OR. (NUMBER ODE 0p eo etic 
LARRAYCI) = -DINF 
GOTO 50 
ELSE 
LARRAYCI) = NUMER/DENOM 
END IF 
IF CLARRAYCI) .LT. 0.0D0) LARRAYCI) = -LARRAYCT) 
CONTINUE 


(Beol - =DinEr 
DO 75 T=1,NCOMP 
IF CLARRAYCI) .Gi. UBESTOaeRHEN 
LBEST = LARRAYCT) 
INDEX = I 
END IF 
CONTINUE 
DE CLBEST .LI = C-CRKVALI DO Sco ORG 
WT) = WIl - WCINDEX, BSPARECINDEX)) + WCINDEX, CBSPARECIHDEX)+1)) 
BOPARECINDEX) = BSPARECINDEX) + 1 
TReeCHII .LE. MAXI) Gomez 
ZHLOW = 0.0D0 
lit = 0.000 
DO 125 T=1,NCOMP 
ZHLOW = ZHLOW + PCI, BSPAREC JD 
WT = WT + WCI,BSPARECI)) 
CONTINUE 


RETURN 
END 
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APPENDIX B 


SAMPLE INPUT FILES 


PRINGEPAL END ITEM FILE 


M6 0A] 5 * 60.0 

EV Te? AL Te ia 40.0 

MIOLAL 5 4 2000.0 
TAS e Es 5a 


PJCC LEAL ENB SLI BM FrbDE SPECIFICATIONS 


LINE COLUMN DATA DESCRIPTION DATA TYPE 
All 7 Nomenclature of end item Character 
S=11 Blank 
12-14 Quantity of end item Integer 
S26 Blank 
19-21 Number of end items 
required operational Integer 
22a Blank 
26-31 Mission duration Real*8 


li 


REPRODUCED AT GOVERNME(Iét FMPRUISE 


SECONDARY REPARABLE COMPONENT FILE 


2815001255367 
2520002241867 
28ToUUZS eco 
2520001549652 
2550000886657 
2520001184942 
292001 UisSs02 
1010010703303 
LOZS5019338S106 
10100 P07 20577 
LO02Z501 07063935 
2615004303480 
252000397555 
2520001943385 
5805014593410 
2520008949552 
2910011719636 
253001 TSG e7 
5865012207848 
4520012055660 
29Z0UUZSI7 276 
2530000886650 
100501 Fees 059 
1005011086434 
1005011457709 


WNW CIR DD PAO AD PO 1 PO PAOD PO PO Re Re RR eRe ee 


20012 
Zuo 
207s. 
CUE 
cone. 
PANE 
alle ake 
20. 
20S. 
20gt, 
2091. 
I25cn 
wee - 
329638 
BAS (e- 
Sao. 
53256 « 
a2 aioe 
a9) 
5258. 
JcIo% 
S2pee 
2007 Ge 
20926. 
2091967. 


on 5 on Ei on 3 em Ei on Fi oo Fe Eom 5 00 Fh 00 8 ee 5 EO EO ROBO RS ROR SE OE OROTOE 


: be f\) 
NR NERA) DR WNO DW OOF NUR ee WN ONS 


_— 
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ZOO 
1000. 
200. 
500. 
10.03 
Z2o0r 
500. 
L000". 
PoO. 
2000. 
Pooo; 
SOG 
OH CHO) 
2000 . 
Zoe 
500. 
1:00: 
PUOO 
100. 
500: 
5/01 0)7- 
2208 
1000. 
2508 
> 00K 


Oooo dao0nd0on0T00onaao nano oo0ooaoo0 0000 0 


_— 


bt 


ee ee 
niilinniooooodMNlilnniininitno OWMNOoOMUuilnilh 


bh he ORO ORR RR eH Orr O 


. 34000 
.00000 
. 00000 
-10110 
.00000 
.00000 
= (0/0) 000 
. 00000 
.00000 
.00000 
ma ehehenone 
»44530 
ste) 3 100) 
.00000 
-20o10 
.00000 
.00000 
.00000 
.0G000 
.00000 
.009000 
.00000 
.00000 
me 100/.01.0 
.00000 


— 
ee ee ee ed ee ee ne.) 


eed mel anal ane el eed ee oe eee) 


6908 
. 0000 
.0000 
. 16am 
0000 
. O00 
.0000 
.0000 
.0000 
.0000 
. 0000 
<.9404 
.0598 
.0000 
. 7348 
. 00 T 
0000 
.0000 
.0000 
.0000 
. OC 0it 
.0000 
OU 
.0000 
.0000 


ee e. 


ee A 


COLUMN 


al 
Pa 1) 


IL 1% 


ISL, 


20-28 


AG 86, 


31-34 


GIES 
a - a2 
43-44 


45-47 


48-49 


30 = 56 


J /—- Je 


a Ob 


TABLE B-2 


DATA DESCRIPTION 





National Stock Number (NSN) 
Blank 

Identification Number 

(li for wisest, 
3 for M1O1A1) 
Blank 


Total exposure time for 
component 


Blank 


Number of observed failures 
for component 


Blank 
Component shipping weight 
Blank 


Maximum possible stockage 
for component 


Blank 

shape parameter of prior 
GEStTrI DOUG On fer seis ie 
failure rate 

Blank 

scale parameter of prior 


distribution for component 
failure rate 
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2 for LVTP7AI1, 


SECONDARY REPARABLE/COMPONENT FILE SPECIFICATIONS 


ye eee i 


Character 


Integer 


eves 


Integer 


Real*8 


Integer 


Real*8 


Real*8 


APRS Neem. 


SAMPLE OUTPUT FILE 
OBJECTIVE FUNCTION VALUE ES ¥=05e820ZGSTso0r00 
THIS SOLUTION IS WITHIN 0.4761377367D+00 % OF THE OPTIMAL 


PCMESSLON SUCCESS e— 30.4139 


REPRODUCED AT GOVERNKIFNI FUE NSE 


NSN END ITEM SPARE LEVEL 
2815001245387 M60Al 2 
2520092241867 M60Al 3 
2815002395819 M60Al 2 
2520001549632 M60Al 2 
2530000886657 M60Al 3 
2520001184942 M60A1 2 
2920010135802 M6 0Al 1 
1010010703803 M60Al 1 
1025019322105 M60Al 5 
1010010720397 M60A1 1 
1025010708993 M60Al l 
2815094303480 LVTP7Al 4 
2520093973334 LVTP7A1 5 
2520001443385 LVTP7AL 4 
5805014593410 LVTP7A1 5 
2520008949532 LVTP7A1 2 
29100117146 36 LVTP7A1 3 
2530011509757 LVTP7A1 2 
5865012207348 LVTP7A1 3 
4520012035660 LVTP7AI 2 
2920002317276 LVTP7A1 1 
2530000886650 LVTP7A1 2 
1005011855059 M101Al 4 
1005011086434 M101Al 5 
1005011457709 M101Al 5 
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APPENDIX D 


GAMS SOURCE CODE 


SETS 
I end item 4/E1875,E03896, E06490/ 
J component 7 1X 257 
K spare stockage level /S0xXS10/ 
L dummy index 7 Det 7; 
SCALARS 


WICAP total weight capacity Cibs) 74464600/ 
PSUC probability of mission completion; 


PARAMETERS 
EQDENSCI) 


Gm? 5 


E0846 
E0640 


PARAMETERS 
NOPERCTI) 
Wee 7 SD 
E0846 
E0640 


PARAMETERS 
DURMISCTI) 
ye 67 5 
E0846 
E0640 


PARAMETERS 
Si ees) 
ZA 


bd bat bt tt bt OO MI ON UN DN NAD 


p~s 
Oo 


iy, 


~<a 
Oo 


density of end item 1 
5 


ee 
37 & 


number end item 1 required operational during mission 


11 
G/ ; 


mission duration for end item 1 Chours rounds miles) 
60 


G0 
20007 ; 


shipping weight of component 53 (lbs) 
2000 
1000 
500 
500 
100 
25.0 
500 
1000 
100 
2000 
1000 
500 
1000 
2000 
250 
500 
100 
1000 
100 
500 
500 
250 
1000 
250 
5007; 


US 


REPRODUCED AT GOVERNAIENT FOCPERIGE 


PARAMETERS 
SPRMAXCJ) 
7 


PO RO ARQ RD A Re OOOH DG 
MSH OwownA i DWIVe © 


MN RO 
in DS Ga 


PARAMETERS 
NOBSCJ) 
wk 


Mm er em OOOO DWP 


ONAN UDWNN HO 


max number of spares allowed for component j 


—— 
min ooonoo MUM OOMoMMNnh 


b— $4 bh 


/3 


number of data points observed for component j 


— 1) 


p— 
“IRN RAR DROWN OO DW OOH WUD Be ee es 0 
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PARAMETERS ; 
ALPHACJ) shape paramtr for Gamma prior distribution 


7a: 0.34 
Z l 
S Ih 
G OMG 
5 Ih 
6 1 
7 ut 
8 i 
9 i 
118 i 
ial I 
2 0.4433 
iL 3 0.6970 
14 ut 
Ss a4) AS 
16 1 
I 1 
18 l 
19 i 
20 it 
rat l 
Ze l 
HES Ih 
24 l 
25 | ae 
PARAMETERS 
BETAC J) scale paramtr for Gamma prior distribution 
Val 58.6903 
2 l 
S if 
G eles 0 
5 it 
6 Ih 
7 jh 
8 i 
9 l 
10 i 
1 i 
V2 142.9404 
is 69.0598 
14 l 
75 214 SoS 
16 IL 
7 1 
18 1 
19 1 
20 l 
Zil i 
Ze 1 
aS i) 
24 i 
25 lea. 
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NS 


REPRODUCED AT GOVERNNEWI!I F 


PARAMETERS 


OKCI,J) equals ] if component 3 iS a component of end item i 
7ZEL67 5 Cle 
F0866 .C1 2X2 ] 
E0640.¢€23x25) Ag 


PARAMETERS 
TXPOSECI) total observed exposure for end item 1 (hrs rds mi) 
PAs eles 165, 209] 
E0846 5256 


E0640 c0I967; 


PARAMETER 


® 
% 


x 
x 


LMTCI,J) max number of failures for component j; 
LMTCI,JI$C OKCI,J) EQ £ ) = EQDENBS(I)D - NGpeas 


Sees 


OKLVLICI,J,K) feasible spare stockage levels for component 3; 
OKLVLICI,J,KI$SCCORDCK)J-1 LE SPRMAXCJ)) AND COKC],J) EQ 1)))) =a 


The following set is necessary to obtain the cumulative distribution 


function of the negative binomial density. 
SES 
ONEVEACl, Jaks eo 
COKLVL2CI,JU,K,LI$C CORDCK)-1 LE SPRMAXCJ)) AND (OKC Eo ee 
CORDCEJ=$1 LE LMICI,JITORDC KD) eKe 


Bro Lae OKLVLZ; 
OPTIONS DECIMALS = 8; 


PARAMETER 
OA, negative binomial paramtr 
QOC1, J) one minus the pe paramtr 
Priel) Pp paramtr raised to power of nobs plus alpha 
THTICI,J,K) total snipping wt of k spares of component j 
Cee, cumulative density for negative binomial 
el hele natural logarithm of cumulative density 
PTMPCI,J,K,L) temporary paramtr needed to sum cdf recursively; 
PCI, J0>SCOKC! Soe?) 


(TXPOSE(I) + BETACJ))/( CEQDENS< I) XDURMIS(I)) + TXPOSECI) + BETACJ)); 
Q(I,JISCOX(1, J) EQ 1) = 1 = PCL, 

PP(I,J)$COKCI,J) EQ 1) = PCI,J) *% CNOBS(J) + ALPHACJ)); 
THTCI,J,K$OKLVLICI,J,K) = SHIPWT(J) * CORDCK) - 1); 
PTMPCI,J,K,"DO")$OKLVLICI,J,K) = 1; 


The following LOOP function is used to calculate the cdf from the 
negative binemial pdf recursively. 
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aoe Cee sede tty doK, L412 )SOKLVL2CI1,J,K,L) = 
een Canosa I tTALP HAC JItTORDCE)—1J/(CORDCL)))XPIMPCIJ,K,L) D3 
Ole s KV SURLVELCI,JS, KD. = 
SUMCLSCCORDCLI~1) LE CLMTCI,JI4+0ORDCK)I-1)), PTMPCI,J,.K,L))*PPCI,J); 
meee OKEVEWCL, I,K) = LOGC CCL,J,K) ); 
DISPLAY F; 
VARIABLES 
ACloyw> KD 1 if stock kK spares for component._3 of end item i 
nee natural log of probability of mission completion; 
BINARY VARIABLE. x; 
EQUATIONS 
WICONST weight constraint equation 
Seeeeiut, J) select exactly [@emfethe possible stockage levels 
MOE measure of effectiveness definition; 


WICONST.. SUMCCI,J,KISOKLVLICI,J,K), THTCI,J,K)¥XCI,J,K)) =L= WTCAP; 


meimeenCieo >< OKCE,J) EQ 1 ).. 
SUM GmEVilCl,J7K), ACI,J,4) ) =E= 1; 


eee. ove (lL, J>RISOKEVEICI,J,K), FCL,J,K)*XCI,J,K) 3 =E= MOE; 
MODEL STOCK /ALL/Z; 

Or rlOns OPTCR = 0.0015, OPTCA = 0; 

SObve otOCK USING MIP MAXIMIZING MOE; 

Bowe = EXPC MOE.L >? 
BoP LAV, Mees, PSUC, 
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